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The extragalactic background light (EBL) is of fundamental importance both for understanding 
the entire process of galaxy evolution and for 7-ray astronomy. However, the overall spectrum 



^ of the EBL between 0. 1 and 1000 jxm has never been determined directly neither from observed 

luminosity functions (LPs), over a wide redshift range, nor from any multiwavelength observa- 
tion of galaxy spectral energy distributions (SEDs). The evolving, overall spectrum of the EBL 
is derived here utilizing a novel method based on observations only. The changing fractions of 



galaxies from redshift 0.2 to 1 are estimated, and two alternative extrapolations of SED types to 
higher redshifts are considered. This allows calculation of the evolving EBL. The EBL uncer- 



> 

quiescent galaxies, star-forming galaxies, starburst galaxies and active galactic nucleus (AGN) 

tainties in our modelling based directly on the data are quantified, and their consequences for 
attenuation of very-high-energy 7-rays due to pair production on the EBL are discussed. It is 
^ concluded that the EBL seems well constrained from the UV to the mid-IR at an intensity level 

roughly matching galaxy count data. Independent efforts from IR and y-ray astronomy are needed 
in order to reduce the uncertainties in the far-IR. 
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1. Introduction 

The formation and evolution of galaxies in the universe are accompanied unavoidably by the 
emission of radiation. All this radiated energy is still streaming through the universe, although 
much is now at longer wavelengths due to redshifting and absorption/re-emission by dust. The pho- 
tons mostly lie in the range of ~0. 1-1000 /im, i.e., ultraviolet (UV), optical and infrared (IR), and 
produce the second-most energetic diffuse background after the cosmic microwave background, 
thus being essential for understanding the full energy balance of the universe. We will account in 
this work for the radiation accumulated by star formation processes through most of the life of the 
universe, plus a contribution from active galactic nuclei (AGNs) to this wavelength range, known 
as the diffuse extragalactic background light (EBL). The direct measurement of the EBL is a very 
difficult task subject to high uncertainties. This is mainly due to the contribution of zodiacal light, 
some orders of magnitude larger than the EBL (e.g., [1, 2]). Other observational approaches set 
reliable lower limits on the EBL, such as measuring the integrated light from discrete extragalactic 
sources (e.g., [3, 4, 5]). On the other hand, there are phenomenological approaches in the literature 
that predict an overall EBL model (i.e., between 0.1 and 1000 /im and for any redshift). These are 
basically of the following four kinds. 

(i) Forward evolution, which begins with cosmological initial conditions and follows a forward 
evolution with time by means of semi-analytical models (SAMs) of galaxy formation (e.g., [6, 
7]. These two papers describe the same modelling but discussed in different contexts). 

(ii) Backward evolution, which begins with existing galaxy populations and extrapolates them 
backwards in time (e.g., [8, 9]). 

(iii) Evolution of the galaxy populations that is infeiTcd over a range of redshifts. The galaxy 
evolution is infeiTcd here using some quantity derived from observations such as the stai^ 
formation rate (SFR) density of the universe (e.g., [10, 11]). 

(iv) Evolution of the galaxy populations that is directly observed over the range of redshifts that 
contribute significantly to the EBL. The methodology presented in Ref. [12], which we term 
empirical, is the first that belongs to this category. The main results presented in that paper 
are summarized here. 

Our EBL estimates from type (iv) methodology are compared with the results from the ap- 
proach type (i) presented in Ref. [6, 7], and with the results from type (ii) methodology by Ref. [9]. 
We consider these two models the most robust among the literature. First, we consider the the- 
oretical approach taken in Refs. [6, 7] as complementary to our observationally motivated one 
to eventually reach a complete understanding of galaxy evolution. Second, approaches type (ii) 
are potentially problematic because they imply extrapolations backwards in time of local or low- 
redshift luminosity functions (LFs). Intrinsically different galaxy populations exist at high red- 
shifts, which cannot be accounted for by these extrapolations. Particularly, [9] use observed LFs in 
the near-IR from the local universe to z = 1.4 for describing the elliptical and spiral populations, 
and only local for describing irregular/starbursting galaxies. They distinguish between these galaxy 
morphologies using images from different satellites. Different local LFs and data sets in the IR are 
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used to constrain the mid and far-IR background. Their modelling is complex and not reproducible. 
Despite these particular problems, this methodology is based upon LFs, quantity directly observed 
and well understood unlike type (iii) models based on parametrizations of the history of the SFR 
density of the universe, quantity with large uncertainties and biases. 

One important application of the EBL for /-ray astronomy is to recover the unattenuated spec- 
tra of extragalactic sources. Our goal is to measure the EBL with enough precision that the un- 
certainties due to the EBL modelling, in these recovered unattenuated spectra, are small compared 
with other effects such as statistical or systematic uncertainties in the 7-ray observations. Examples 
of this are discussed in Sec. 4. 

Throughout this work, a standard ACDM cosmology is assumed, with matter density Q,,n = 
0.3, vacuum energy density Q.^ = 0.7, and Hubble constant Hq = 70 km s^'Mpc^^. 

2. Methodology 

Our model is based on the rest-frame K-band galaxy LF in Ref. [13] and on multiwavelength 
galaxy data from the All-wavelength Extended Groth Strip International Survey (AEGIS \ [14]) of 
about 6000 galaxies in the redshift range of 0.2-1. These data sets are put together in a very trans- 
parent and consistent framework. The [13] LF is used to count galaxies (and therefore to normalize 
the total EBL spectral intensity) at each redshift. The LF as well as our galaxy sample is divided 
into three magnitude bins according to the absolute rest-frame K-band magnitude, i.e., faint, mid- 
dle and bright. Within every magnitude bin, an SED type is statistically attached to each galaxy 
in the LF assuming SED-type fractions that are a function of redshift within those magnitude bins. 
This is estimated by fitting our AEGIS galaxy sample to the 25 galaxy-SED templates from the 
SWIRE^ library [15]. Then, luminosity densities are calculated from these magnitude bins from 
every galaxy population at all wavelengths, and finally all the light at all redshifts is added up to 
get the overall EBL spectrum. In this approach, it is possible to directly calculate the contribution 
to the EBL from all redshift bins, as well as the evolution of the EBL spectrum with redshift and 
the processes related to this evolution, by sources of all the 25 SED types considered. 

A multiwavelength galaxy catalogue built from AEGIS for this work is used. This catalogue 
contains 5986 galaxies, all in the Extended Groth Strip (EGS). It is required that every galaxy in the 
sample have 5a detections in the B, R, I, Ks and Infrared Array Camera (IRAC) 1 bands, and ob- 
servations (but not necessarily detections) in the IRAC 2, 3, 4 and Multiband Imaging Photometer 
for Spitzer (MIPS) 24 bands. These 5a upper limits are given by the following fluxes: 1.2, 6.3, 6.9 
and 30 /iJy for IRAC 2, 3, 4 and MIPS 24, respectively, according to Ref. [16] for the IRAC bands 
and Ref. [17] for MIPS 24. In addition, 1129 of these galaxies have Galaxy Evolution Explorer 
(GALEX) detections in the far-UV and 2345 galaxies in the near-UV. In our sample, 4376 galaxies 
have the highest quality spectroscopic redshifts measured by the Deep Evolutionary Exploratory 
Probe 2 team (DEEP2 DR3, [18]), with the Deep Imaging Multi-Object Spectrograph (DEIMOS) 
spectrograph [19] on the Keck II telescope in an area of about 0.7 deg^ in the sky. All the other 
galaxies in the sample (1610 galaxies) have secure photometric redshifts, more than 80 per cent 
with uncertainty in redshift less than 0. 1 . 

' http://aegis.ucolick.org/ 

^http://www.iasf-milano.inaf.it/'~polletta/templates/swire_templates.html 



3 



EBL inferred from AEGIS galaxy SED-type fractions 



A. Dommguez 



X Slat-lot mi tig 




redshift 

Figure 1: Galaxy SED-type fractions from our catalogue (after the Xred cuts) of the different populations 
versus redshift according to our multiwavelength fits. We mark with crosses our fractions from z = 0.9 — 0.3. 
The circles at z = 0.1 are fractions computed from the SDSS-based sample (see Ref. [12] for details). We 
show with a shadow area the uncertainties from our lower limit for the errors as well as for our x^ed '^^^ for 
fits. The uncertainties are around ±0.1. 



3. Results 

3.1 Galaxy SED-type fractions 

The Le PHARE code is used to fit every galaxy in our sample to the 25 SWIRE templates. 
For clarity, we will compress in our discussion (but not in our calculations, where they will remain 
independent) the 25 SED types in the SWIRE library to four groups: quiescent galaxies, star- 
forming galaxies, starbursts and AGN galaxies. 

To avoid accounting for bad fits, which do not correctly describe the galaxy photometric data, 
a cut in Xred ~ is applied, with given by Le PHARE and n degrees of freedom (bands 
with detections). We have checked carefully that x}ed < 30 is a red good value for quiescent, star- 
forming and starburst galaxies, but AGN galaxies are systematically worse fits, probably due to the 
fact that there is a large range in AGN SED shapes due to multiple emission components which 
cannot be easily encapsulated in a few templates and a stronger cut is needed xled — 

Fig. 1 shows the galaxy-SED-type fractions for four different redshift bins up to z = 1, where 
we have chosen bins of z = 0.2 for statistical reasons. The redshift range shown corresponds almost 
to 60 per cent of the age of the universe. The shadow regions are the uncertainties due to the lower 
limits on the photometric errors for the catalogue and for the Xred "^^^^ region is calculated 
changing the lower limits from 1 to 10 per cent in steps of 1 per cent and applying extreme cases for 
the cuts for every lower limit. The boundaries from these calculations lead to the shadow regions. 
The fractions adopted for the model are marked with crosses and wider lines. We observe that the 
fraction of quiescent galaxies (dashed-red line) increases by a factor of ~ 2 from z ~ 0.9 to 0. 3, 
while the star-forming fraction (dotted-dashed-blue line) keeps roughly constant for the full redshift 
range peaking at z = 0.5. Starburst-type galaxies (dotted-green line) decrease very quickly from 
z ~ 0.9 and reach almost at z ~ 0.5. On the other hand, the AGN-type fraction (solid-gray line) 
is roughly constant from z ~ 0.9 to 0.7 and then decreases to 0.02 at z ~ 0.3. This result should not 
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be considered a complete picture of the evolution of the galaxy populations in the universe since 
these fractions depend on the color-magnitude limits of the survey. But what is certainly described 
is the population of galaxies that contribute the most to the EBL around the knee of the LF. 

For the high-redshift universe (z > 1, where there are no galaxies in our sample), two different 
cases are considered for the evolution of the galaxy-SED-type fraction. It is shown that our results 
are not changed significantly except in the far-IR by these two choices. For the redshifts less than 
those of the most distant known 7-ray sources, and redshifts where sources are likely to be found 
in the near future by Imaging Atmospheric Cherenkov Telescopes (lACTs), we find that there is al- 
most no change in the EBL even with a fairly large adjustment in the evolution of galaxy-SED-type 
fractions. The fiducial choice is to keep constant the fractions computed for our highest redshift 
bin. This choice is made for simplicity, due to the difficulty in the multiwavelength classification 
of distant galaxies with current instruments. As an alternative approach, we choose to increase 
hnearly with redshift the starburst-hke fraction from our calculated 16 per cent at z = 0.9 up to 60 
per cent at z = 2, while decreasing at the same rate the quiescent and star-forming galaxies. The 
weight of each of the 25 SWIRE templates is changed in the same proportion. The fractions are 
kept constant at z = 2 for z > 2. This approach is called high starburst and it is used to determine a 
likely upper limit on the EBL at long wavelengths 

3.2 Extragalactic background light 

The local galaxy luminosity density, its evolution over redshift at different wavelengths for the 
two extrapolations previously considered for the high-redshift fractions, as well as an estimation 
of the SFR density of the universe are studied in Ref. [12]. Other quantities such as the EBL 
evolution^ are discussed there as well. Here, we only show in Fig. 2 the EBL in the local universe, 
with its uncertainties, compared with direct and indirect observational data, and other EBL models. 
Fig. 2 suggests that the EBL coming from galaxies is already well constrained in the region from 
the UV up to the mid-IR but not in the far-IR. Galaxy counts from very deep surveys taken with 
very sensitive instruments [3, 4, 5] should be considered as a good estimation of the true EBL from 
galaxies. On the other hand, different fully independent modelings based on different galaxy data 
sets [7, 9, 12] agree in the specific intensity level of the EBL. In particular, galaxy count data ai"e in 
excellent agreement with our results. 

4. Application to y-ray attenuation 

The EBL has important implications for the interpretation of data taken using recent 7-ray 
experiments such as the Fermi satelhte and lACTs (VERITAS, HESS, and MAGIC) due to the 
photon-photon pair production between 7-ray photons traveling across cosmological distances and 
EBL photons [20, 21]. Blazars are an important source of extragalactic 7-ray emission and have 
become a relevant tool for indirectly measuring the EBL. These objects are believed to be an ex- 
treme category of AGNs. Their emission, which occurs at all wavelengths of the electromagnetic 
spectrum, comes from super-massive black holes (with masses of > 10^ M0) swallowing matter 
accreted from their surroundings. In general, AGNs are characterized by a beamed emission per- 
pendicular to the accretion disc known as jets, which are pointing towards us in the case of blazars. 

^EBL specific intensities are publicly available at http://side.iaa.es/EBL 
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Figure 2: The solid-black line is the extragalactic background light calculated by the fiducial extrapolation 
of the galaxy SED-type fractions for z > 1 . Uncertainties in the our EBL estimation are shown with a shadow 
area (see Ref. [12] for a discussion on this and for details on the references). The envelope of the shadow 
region within the dashed line at wavelengths above 24 /im shows the region where there is no photometry in 
our galaxy catalogue. 



Following both theoretical arguments [22, 23] and observational facts [24, 25], it is assumed 
that no intrinsic (or EBL-coiTccted) VHE spectra from blazars might be fitted to a power-law with 
indexes hai^der than 1.5. We now proceed to test in Fig. 3 whether the observed spectra of three new 
measurements of high-redshift AGNs (other different spectra were considered in Ref. [12]) satisfy 
the condition that the intrinsic spectrum corrected by the attenuation derived with our EBL model'* 
has Tint > 1.5. We consider the following blazars: 3C 66A at z = 0.444 observed by MAGIC 
[26], 3C 279 at z = 0.536 observed by MAGIC in the 2007 observational campaign [27], and the 
discovery of PKS 1222-1-216 in the VHE regime [28], the second most distant flat-spectrum radio 
quasar known (z = 0.432). These three blazars are plotted in Fig. 3, where the legends show that 
the condition F,,,, > 1.5 is satisfied. We note that in the 3C 279 case, only having three data points 
makes the fit no statistically reliable. 

It is confirmed from the study of these blazai^s the conclusions obtained in Ref. [12]. First, our 
EBL is generally compatible with the expected hardness of the EBL-corrected slopes. However, it 
is clear that a simple SSC model cannot explain any flatness at the highest energies of the EBL- 
corrected spectra of 3C 66A, which suggests that some extension to the model may be necessary 
such as an external photon region, a better understanding of the lACT systematic uncertainties or 
even a revision of the propagation mechanisms mainly through the intergalactic medium [29]. 

Second, the uncertainties in the EBL-corrected spectra are dominated by other effects different 
than EBL modelling as shown in the index uncertainties in Fig. 3. 



Optical depths are publicly available at http://side.iaa.es/EBL 



6 



EBL inferred from AEGIS galaxy SED-type fractions 



A. Dommguez 




Figure 3: VHE spectra observed (blue) and EBL corrected (red) from the attenuation calculated with our 
EBL model of three high-redshift blazars measured by MAGIC [26, 27, 28]. The straight-red line is the 
best-fitting power law for every blazar with index r„„. Uncertainties from the EBL modelling as well as 
statistical and systematic errors are shown with a shadow region. The first uncertainties in the index are due 
to the EBL modelling and the second uncertainties are statistical plus systematic errors. Intrinsic energies 
are slightly shifted for clarity. 



5. Conclusions 

A novel, robust, and powerful method based on observations to derive the evolving spectrum 
of the EBL between 0.1 and 1000 /^m was thoroughly presented in Ref. [12] and reviewed here. 
This model is based on the observed rest- frame K-band galaxy LF over redshift found in Ref. [13], 
combined with an estimation of galaxy-SED-type fractions based on a multiwavelength sample of 
~ 6000 galaxies from AEGIS. This model has the following main advantages over other existing 
EBL models: transparent methodology, reproducibility, and utilizing direct galaxy data. The best 
available data sets are used ([13]'s LF and the AEGIS galaxy catalogue) observed over a wide 
redshift range. The galaxy evolution is directly observed in the rest-frame K band up to z = 4. 
Observed galaxies up to z = 1 from the UV up to 24 /xm with SEDs of 25 different types (from 
quiescent to rapidly star-forming galaxies and including AGN galaxies) are taken into account 
in the same observational framework. A study of the uncertainties to the model directly from 
the data (such as uncertainties in the Schechter parameters of the [13] LF and the errors in the 
photometric catalogue) is done, and their propagated uncertainties to the /-ray attenuation are 
studied. Two extrapolations of the galaxy-SED-type fractions to z > 1 were considered, showing 
that these assumptions only affect the far-IR, where the uncertainties of the modelling are the largest 
because of the lack in our catalogue of far-IR photometry and the poor understanding on galaxy 
SED at large redshifts. 

We concluded that the EBL from galaxies seems already well constrained from UV to mid-IR 
wavelengths, even though uncertainties are still large in the far-IR. Furthermore, discoveries of 7- 
ray from distant blazars (e.g., [26, 27, 28]) support the EBL specific intensity level derived from 
galaxy count and recent EBL models [7, 9, 12]. As discussed the VHE recovered spectra of blazars 
up to ~ 10 TeV are dominated by statistical and systematic uncertainties in the observations rather 
than EBL modelling. We highlight that the EBL specific intensity calculated with our method is 
matching the lower limits from galaxy counts, which implies the highest transparency of the uni- 
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verse to y-ray allowed by standard physics. This predicts a promising future for the new generation 
of lACTs, namely CTA. 
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